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We revisit the critical behavior of the sub-ohmic spin-boson model. Analysis of both the leading 
and subleading terms in the temperature dependence of the inverse static local spin susceptibility 
at the quantum critical point, calculated using a numerical renormalization-group method, provides 
evidence that the quantum critical point is interacting in cases where the quantum-to-classical 
mapping would predict mean-field behavior. The subleading term is shown to be consistent with 
an ui/T scaling of the local dynamical susceptibility, as is the leading term. The frequency and 
temperature dependences of the local spin susceptibility in the strong-coupling (delocalized) regime 
are also presented. We attribute the violation of the quantum-to-classical mapping to a Berry-phase 
term in a continuum path-integral representation of the model. This effect connects the behavior 
discussed here with its counterparts in models with continuous spin symmetry. 
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I. INTRODUCTION 

The standard Kondo model, describing exchange scat- 
tering between a local moment and a fermionic bandji 
can be generalized to various Bose-Fcrmi Kondo mod- 
els that include coupling of the impurity to one or more 
bosonic baths. It is conventional to consider bosonic 
baths described by a power-law density of states 
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where w p is the dispersion of the bosonic bath (see be- 
low), is the Heavisidc function, and uj c is a high-energy 
cut-off. For a subset of sub-ohmic (positivc-e) baths 
corresponding to < e < 1, the Bose- Fermi Kondo 
problem has a quantum (temperature T = 0) phase 
transition^— between a Kondo phase in which the ground 
state is a Kondo singlet formed between the local mo- 
ment and conduction-electron spins, and a local-moment 
phase in which the coupling to the bosonic bath inhibits 
the Kondo effect. In all the Bose-Fermi Kondo mod- 
els that have been studied — corresponding to Ising, XY, 
and SU(2) symmetry of the bosonic couplings — this tran- 
sition is continuous. 

The Landau theory for continuous phase transitions 
invokes the fluctuations of an order parameter. For the 
Bose-Fermi Kondo problem, a suitable order parameter 
is the local magnetization M = (S z ), which vanishes in 
the Kondo phase but is nonzero throughout the local- 
moment phase. Within the Landau approach, the quan- 
tum criticality is described by a local 4 theory with a dis- 
sipative quadratic term corresponding to a long-ranged 
interaction in imaginary time: ^S z (t)xo ( r ~ T ')S Z ( T ') 
with 
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In keeping with the standard theory of quantum crit- 
icality, this description is called a quantum-to-classical 
mapping. 

In a study of the Bose-Fermi Kondo model with SU(A?") 
spin symmetry and SU(M) channel symmetry, it was 
shown in the limit of large N and M with M/N fixed that 
the quantum critical point (QCP) breaks the quantum- 
to-classical mapping.— The violation is especially clear 
for i < e < 1. In this range, the local-0 4 description 
gives rise to a Gaussian fixed point r- and a dangerously 
irrelevant </> 4 coupling leads to a critical static spin sus- 
ceptibility 
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However, the large- N results for the local spin suscepti- 
bility and Green's function in the quantum critical regime 
of the SU(iV) x SU(M) Bose-Fermi Kondo model satisfy 
oj/T scaling, which implies that the QCP is interacting. 
In particular, 
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hold over the entire range < e < I . This violation of the 
quantum-to-classical mapping survives the inclusion of 
1/N corrections, and has been attributed to the effect of 
the Berry phase in the spin path-integral representation 
of the local moment 4 

Subsequently, a violation of quantum to classical map- 
ping was also discussed for the sub-ohmic spin-boson 
model.— This model (fully specified in Section II) couples 
the component S z of an SU(2) local spin to the displace- 
ment of a sub-ohmic bosonic bath with coupling constant 
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g and the component S x to a transverse magnetic field 
r. The model has a QCP between a delocalized phase 
in which the impurity degree of freedom is quenched by 
the field and a boson-dominated localized phase that re- 
tains a two-fold local-moment degree of freedom. These 
phases are analogous to the Kondo and local-moment 
phases, respectively, of the Bosc-Fcrmi Kondo model. In 
the special case where the bosons form an ohmic bath, 
the spin-boson model and the fermionic Kondo model can 
be transformed into one another via bosonizationiS; how- 
ever, this transformation breaks down for e > 0. In Ref. 
9, it was shown using a bosonic extensio n 11 ' 12 of the nu- 
merical renormalization-group (NRG) method that the 
critical exponents for | < e < 1 obey hyperscaling (as 
they do for < e < i). In particular, just as for the 
Bose-Fermi Kondo model^i 5 - the critical local susceptibil- 
ity satisfies Eqs. (@| and JS]), suggestive of an interacting 
fixed point with to /T scaling — a conclusion at odds with 
the quantum-to-classical mapping. 

Recently, two sets of Monte-Carlo calculation s) 5 ' 13 were 
carried out for a one-dimensional classical spin chain with 
long-range Ising interactions J2i j Jij^iSj, where Jij ~ 

1/K — r j\ 2 ~ e f° r \ r i ~ r j\ ^> r o- Referenced applied a 
cluster Monte-Carlo method at various nonzero values of 
the short-range cutoff ro, while Ref. UH used a similar 
method after taking the limit To —> 0. These works both 
confirmed an earlier conclusion^ that the critical points 
of such long-ranged classical spin chains are interacting 
for < e < \ and Gaussian for 4 < e < 1, consistent 
with the predictions of the local cf) 4 theory.— ^ 

The results of the classical Monte-Carlo calculations 
have been interpreted in two very different ways. In 
Ref. [fl we suggested that these results differ from NRG 
calculations for the Bose-Fermi Kondo model because a 
classical spin chain of length L sites does not faithfully 
represent the quantum-mechanical model at temperature 
T = 1/(Ltq). More specifically, the limit ro — > of the 
classical spin chain does not reproduce the path integral 
of the quantum problem because this limit smears the 
topological effect of the Kondo spin flips. Along this line, 
it has recently been showni^ that a proper path-integral 
for the sub-ohmic spin-boson model involves a Berry- 
phase term. Numerically, Ref. [f| generalizing the NRG 
result of Ref. 0, showed that the temperature dependence 
of the local spin susceptibility of the Bose-Fermi Kondo 
model obeys Eq. (j4]) over about 20 decades in tempera- 
ture. 

By contrast, Ref. [TH assumed the the classical spin 
chain faithfully represents the spin-boson model. Cor- 
respondingly, it interpreted the Monte-Carlo result 
as demonstrating a fundamental error in the NRG 
results ^ 11 ' 12 for the sub-ohmic spin-boson model. It was 
further suggested in Refs. [T3 and [H that over the range 
5 < e < f, the NRG result for the temperature depen- 
dence of the local spin susceptibility obeys Eq. © due to 
an artifact of the method, reminiscent of the effect of a 
temperature-dependent mass terrn^ discussed in Ref. [H, 
and that removal of a spurious T 1 ~ e term from x~j t (T) 



exposes an underlying T 1 / 2 term representing the true 
critical behavior. 

In this paper, we investigate more thoroughly the tem- 
perature and frequency dependences of the local spin sus- 
ceptibility of the spin-boson model, as calculated using 
the bosonic NRG approach i 11 ' 12 The subleading term in 
the temperature dependence at lu — is shown to be de- 
scribed by an exponent that depends on e and exceeds 
h , contradicting a central assumption of Ref. [IH The ex- 
ponent 2/2 of the subleading term in the frequency depen- 
dence at T = satisfies yi = X2, paralleling the equality 
y = x of the leading exponents. This provides evidence 
that both the leading and subleading terms in the lo- 
cal susceptibility satisfy ujJT scaling. Finally, we show 
that NRG gives consistent frequency and temperature 
dependences of the local spin susceptibility in the strong- 
coupling (delocalized) phase. All these features indicate 
that the temperature and frequency dependences of the 
critical local spin susceptibility are consistent with an 
interacting fixed point for | < e < 1, thereby augment- 
ing previous evidence for the violation of quantum-to- 
classical mapping in the sub-ohmic spin-boson and Bose- 
Fermi Kondo models. 

The remainder of the paper is organized as follows. 
The model and the bosonic NRG method are described 
in Sec. HI1 Section Mil is devoted to the analysis of the 
leading and subleading terms in the temperature and fre- 
quency dependences of the inverse local spin susceptibil- 
ity at the QCP. In Sec. IIV1 we discuss the local suscepti- 
bility in the strong-coupling (delocalized) phase. Section 
fVl addresses the implications of our results, the form of 
the proper path integral for the spin-boson model, and 
the role of the Berry phase. The paper concludes with a 
brief summary in Sec. EH 



II. MODEL AND SOLUTION METHOD 

The spin-boson model is described by the Hamiltonian 

H SB = -TS X + gS z J2 (Vp + <t>± P ) + X>*> <^ > 

V P 

where S a = er Q /2 with a a being a Pauli matrix, and <fi p 
annihilates a boson in a bath whose density of states is 
specified by Eq. ([1]). The cutoff frequency uj c entering 
Eq. (p} will be set to unity and will henceforth serve as 
the unit of energy. 

We have studied the model using the NRG method 
described in Refs. [ll| and [TH- The bosonic bath is di- 
vided into a set of bins spanning oscillator frequencies 
A-( fe +!) < u) < A- k for A > I and k = I, 2, 3, 
Within each bin, the continuum of bath states is replaced 
by a single state, namely, the linear combination of states 
that couples to the impurity spin. Then the Lanczos 
method is used to map the spin-boson Hamiltonian to 
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H n = -TS x + gS z J2(b + b > 



^2 \e m b\ n b m + t m {b\ n b m _ 1 



H.c. 



(7) 



As a result of the logarithmic binning, the tight-binding 
coefficients e m and t m that encode the bath density of 
states decay as A~ m . This decay allows the Hamilto- 
nian H^ G to be solved iteratively starting with H and 
using the eigensolution of H n to construct the basis of 
H n+ \. The basis of each site of the bosonic chain must 
be truncated at b^b^ < N h . Even with restriction, the 
dimension of the Fock space grows exponentially with it- 
eration number n. After a few iterations, it is possible 
to retain only the N s many-body eigenstates of lowest 
energy after iteration n. All the results reported below 
were obtained using Nb = 24 and N s = 300. 

Our focus in this paper is on the temperature and 
frequency dependences of the local spin susceptibility. 
Within the NRG approach, the local static susceptibility 
is calculated as 
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where ft. is a local magnetic field coupling to the S z com- 
ponent of the localized spin through an additional Hamil- 
tonian term AH = hS z . The imaginary part of the local 
dynamical susceptibility can be computed as 



x 5(uj - E nM +E njj ), 
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where \j) n is a many-body eigenstatc of iteration n with 
energy E n j, and Z n = J2j e~ E "-^ T . In order to mini- 
mize known artifacts of the NRG discretization, we cal- 
culate the local susceptibility only at iterations of the 
same parity (either n even or n odd) and at only one 
temperature or frequency per iteration: T n ,oj n oc A~™. 



III. LOCAL CORRELATION FUNCTIONS AT 
THE QUANTUM CRITICAL POINT 

The QCP is located in our NRG calculations by fixing 
the transverse field T and tuning g. The critical strength 
g c is identified as the value of g at which the scaled NRG 
many-body energies A n E n are independent of the itera- 
tion number n, signifying that the system is at a scale- 
invariant fixed point. This invariance is illustrated in 
Fig. [I] for two values of e. The critical many-body spec- 
trum is qualitatively similar for all baths corresponding 
to < e < 1. 




FIG. 1: (Color online) NRG spectra of the spin-boson model 
at its QCP: Scaled NRG energy A n E n vs iteration number 
n for a transverse field V = 0.01, NRG discretization A = 3, 
and bath exponent (a) e = 0.8, (b) e = 0.2. The flatness in 
the n-dependence of the scaled energy shows that the system 
is at a renormalization-group fixed point; here, it corresponds 
to the QCP at g = g c . 



A. Temperature dependence of the critical local 
susceptibility 

At the QCP, g = g c , the local static spin susceptibility 
has a singular temperature dependence. Figure [Ha) is a 
log-log plot of vs T for e = 0.8. A least-squares fit 
(dashed line) of the logarithm of x _1 to the logarithm of 
a simple power law 



f(T) = aT* 
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over 10 18 < T < 10 yields a critical exponent of 
x = 0.202. Similar behavior is observed in the local static 
susceptibility for e = 0.7, fitted over 10~ 19 < T < 10~ 4 
by x = 0.301 [Fig. Eta)]; for e = 0.6, fitted over 10~ 17 < 
T < 10~ 4 by x = 0.406 [Fig. Ufa)]; and for e = 0.2, fitted 
over 10~ 13 < T < 10~ 4 by x = 0.809 [Fig.EJa)]. These 
temperature exponents are consistent with x = 1 — e to 
within about 1.5%. 

The modest discrepancies between the exponents x ex- 
tracted from the numerical data and their interacting val- 
ues 1 — e can be attributed to a combination of fitting un- 
certainties and errors associated with the NRG method. 
Any NRG calculation contains discretization errors in- 
troduced by working with A > 1, as well as truncation 
errors arising from retaining only N s many-body states 
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FIG. 2: (Color online) Temperature dependence of the local 
static susceptibility at the QCP (g = <? c ). (a) x' 1 ^) for 
e = 0.8, r = 0.01, and A = 3. The dashed line represents a 
fit of x _1 ( T ) to tne leading term f(T) specified in Eq. (fTOl) . 
which yields an exponent x = 0.202. The solid line is a fit 
to Eq. in terms of both leading and subleading terms, 
(b) Residual x _1 ( r ) - f( T ) after subtraction of the fitted 
leading term from (a). The fitted solid line gives a subleading 
exponent X2 = 0.852. The residual clearly does not have a 
T 1//2 dependence (dashed line). 

after each iteration. For the bosonic bath treated here, 
there is also the need to truncate the dimension of the 
Fock space of each bosonic orbital on the NRG chain 
to a finite value A^. One should always bear in mind 
such sources of systematic error, as one would for any 
numerical method. Still, several features of the leading- 
order results serve as nontrivial checks; in particular, the 
NRG value of the frequency exponent y of the critical 
local susceptibility (see Sec. IIII B[) agrees with that of an 
e expansion^, and (b) the temperature exponent x and 
several other critical exponents satisfy hypcrscaling rela- 
tions (as discussed in Sec. |V|) . 

Motivated by considerations laid out in the introduc- 
tion, we further explore the critical behavior by exam- 
ining the subleading temperature dependence defined 
through the fitting function 

X -\T) = aT x + a 2 T x \ (11) 

For the case e = 0.8 shown in Fig.^a), fitting to Eq. ([IT]) 
yields a subleading exponent x 2 = 0.852. The quality of 
the fit of the subleading term can be seen more clearly in 
Fig. CUb), which plots the residual x" 1 ^) - f(T) after 
subtraction of the fitted leading term. The residual is 



FIG. 3: (Color online) Like Fig.rjbut for e = 0.7, V = 1, and 
A = 4. 

seen to be well-described by a power law over about four 
decades of temperature 10~ 6 < T < 10~ 2 . The fitted 
exponent is different from the value x 2 — \ that would 
correspond to the dashed line shown in Fig.[2jb). 

Similar results hold for e = 0.7, fitted over 10~ 5 < 
T < lO^ 1 by x 2 = 0.671 [Fig. Mb)}; for e = 0.6, fitted 
over 10~ 5 < T < 0.5 by x 2 = 0.665 [Fig. gjb)]; and for 
e = 0.2, fitted over 5 x 10~ 5 < T < 0.5 by x 2 = 1.25 [Fig. 
E^b)] . For each of the three cases in the range \ < e < 1 
where the validity of the quantum-to-classical mapping 
is at issue, we find no evidence for the T 1 / 2 contribution 
to that is assumed in Ref. [H. 

B. Frequency dependence of the critical local 
susceptibility 

Next, we investigate the frequency dependence of the 
zero-temperature critical local dynamical susceptibility 
x(w). Since the NRG delivers x"( w ) = I m x( w ): this 
quantity is the focus of our analysis. If necessary, x'^) = 
Rex(w) can be obtained via a Hilbcrt transform of x"(w), 
but this procedure introduces some error due to the lower 
accuracy of the NRG-calculated x"( w ) a t high frequen- 
cies. 

FigurelUJa) plots x"( w ) f° r e = 0-7- The low-frequency 
behavior is singular: 



The exponent is, to an accuracy considerably better than 
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FIG. 4: (Color online) Like Fig. H but for e = 0.6. 



1%, y = 1 — e = 0.3. This implies that the full retarded 
susceptibility has the form specified by Eq. ((5). The 
equality of y and the temperature exponent x is con- 
sistent with an ui/T scaling form for the leading scaling 
term in the inverse local susceptibility. 

We will assume that x(w) has a subleading term with 
exponent y 2 defined through 
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Correspondingly, for uj > 0, 
1 
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where b = B / sin(7ry/2), 

"cos(7r(y + y 2 )/2) , sm(7ry 2 /2) 
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and the approximation in Eq. (1141) involves ignoring 
higher-order terms ~ uj 2y2 ~ v . Fitting to Eq. (fT4| [Fig. 
in^b)] yields an exponent y 2 = 0.73 that agrees with the 
subleading temperature exponent x 2 = 0.67 [Fig. [3jd)] 
to within about 9%. Given that extracting y 2 from 
l/x"(w) is less accurate than determining the exponent 
from l/x(w) [whose determination requires a calculation 
of x'( w )i however], we interpret our results as being con- 
sistent with y 2 = x 2 . This suggests that the subleading 
term of the inverse local spin susceptibility also obeys 
lj/T scaling. 



FIG. 5: (Color online) Like Fig. H but for e = 0.2. 
(a) io 8 




FIG. 6: (Color online) Frequency dependence of the local 
dynamical susceptibility at the QCP (T = 0, g = g c ). (a) 
Inrx(i^) for e = 0.7, F = 1, and A = 4. The solid line repre- 
sents a power-law fit to Eq. (|12[) . which yields the exponent 
y — 0.3002. (b) Residual [Imx(aj)]" 1 — buj v after subtraction 
of the leading term inferred from (a). The fitted solid line 
gives a subleading exponent = 0.729. 
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FIG. 7: (Color online) Like Fig. [6] but for e = 0.2, T = 0.01, 
and A = 3. 
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FIG. 9: (Color online) (a) Inverse X _1 (T) of the local static 
susceptibility at different bosonic couplings g < g c or, equiva- 
lently, 5g = {g c -g)/gc > 0, for (a) e = 0.7, F = 1, and A = 4; 
and (b) e = 0.2, T = 0.01, and A = 3. 
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The same procedure can also be applied to the case 
e = 0.2, as shown in Figs.JTJa) andJTJb). The extracted y 
is again to a good accuracy equal to 1 — e = 0.8, the same 
as the leading temperature exponent x, so the leading 
term of the inverse local susceptibility is consistent with 
uj/T scaling. The fitted value of the subleading frequency 
exponent is yi = 1.32, which agrees with the subleading 
temperature exponent xi to within 6%. It is worth noting 
that for this value of e, which falls in the range < e < |, 
the level of agreement between X2 and y2 is similar in 
percentage terms to that for e = 0.7 in the range | < 
e < 1. 
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FIG. 8: (Color online) Imaginary part Imx(w) of the zero- 
temperature local dynamical susceptibility at four different 
bosonic couplings g < g c or, equivalently, 8g = {gc—g)/gc > 0, 
for (a) e = 0.7, T = 1, and A = 4; and (b) e = 0.2, T = 0.01, 
and A = 3. 



IV. LOCAL CORRELATION FUNCTIONS IN 
THE DELOCALIZED PHASE 

To further analyze the NRG results, we turn to the lo- 
cal susceptibility in the delocalized phase, g < <? c , which 
has not received much attention in previous studies. Fig- 
ure [S^ a) shows the u> dependence of x"( w ) a t zero tem- 
perature for e = 0.7 and several values of Sg = (g c ~g) / gc- 
The leading frequency dependence is seen to be x"(w) ~ 
|w| 1_e . The results are well described by 

X(T = 0,u,g < g c ) 



A + B(-iuy 



(16) 



The decrease of x(u; = 0) with increasing 5g (i.e., de- 
creasing g) in Fig. (5Ja) reflects the increase of A as g 
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moves away from g c . These conclusions are also valid for 
e = 0.2, as seen in Fig. |5Jb). 

Figure [9^a) shows that for e = 0.7, the local static 
susceptibility in the delocalized phase is well described 
by 

X (T,u = Q,g<g c ) = A+ * yl _ 6 - (17) 

The same conclusions, once again, apply to e = 0.2, as 
seen in Fig. [9jb). 

We close this section with two general remarks. First, 
the exponent 1 — e describing the frequency dependence 
(and, by extension, the temperature dependence), takes 
the largest possible value that satisfies the Griffiths 
inequality— which states that correlation functions de- 
cay in time no faster than does the interaction [as spec- 
ified in this case by Eqs. (f2"rJ)) and (j2"T|) below]. Second, 
this frequency and temperature dependence only exists in 
the delocalized phase for g =/= 0, implying that turning on 
the coupling g to the bosonic bath within the delocalized 
phase is a singular perturbation. 

V. DISCUSSION 
A. Implications of our results 

We have considered in detail the critical behavior of 
the local spin susceptibility of the sub-ohmic spin-boson 
model calculated using the NRG. For four different bath 
exponents spanning the range 0.2 < e < 0.8, we have 
confirmed that the leading term in the inverse static lo- 
cal spin susceptibility has the T 1_e temperature depen- 
dence described by Eq. (fj| and demonstrated that the 
subleading term varies as T X2 with X2 > \. 

This analysis allows us to assess the alternative inter- 
pretations (as outlined in the Introduction) of the Monte- 
Carlo results for the critical behavior of the local suscepti- 
bility in a long-ranged classical Ising chain. In the regime 
| < e < 1, the chain yields a temperature exponent of 
■i. This originates from a dangerously irrelevant vari- 
able, which also leads to a violation of uj/T scaling. In 
particular, the T 2 dependence of of the inverse local sus- 
ceptibility is not accompanied by a (— iu>)^ dependence 
on frequency. These results differ from those of NRG 
calculations for the quantum model. 

The NRG results presented in Sec. Mil arc consistent 
with the interpretation advanced in Ref. [5| in the context 
of the Bose-Fermi Kondo model, namely the difference 
of the Monte-Carlo results from the NRG results and 
reflects a violation of the quantum-to-classical mapping. 
The latter is the result of a Berry-phase term, which we 
discuss in the next subsection. 

Our results are inconsistent with the interpretation put 
forward in Refs. [TH, [I?] and [HI, which advocates the va- 
lidity of quantum-to-classical mapping. In this picture, 
the intrinsic temperature dependence of the inverse crit- 
ical local susceptibility of the NRG calculation should be 
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FIG. 10: Schematic log-log plot of the inverse static local 
susceptibility x _1 vs temperature as obtained using the NRG 
(cf. Figs. [2H2] of the present paper) and NRG* [cf. Figs. 7(b) 
and 8(b) of Ref.—] procedures. The curvature of the NRG* 
line suggests that a T 1 ' 2 dependence at T < T* crosses over 
at T > T* to a T z dependence with a z smaller than ^; this in 
turn suggests that the asymptotic low-temperature behavior 
has a temperature exponent less than |. 



Ta, which is masked by an artificial leading T 1_e term. 
We have shown here that there is no subleading term 
in the NRG results. 

Reference [H suggested the possibility of removing the 
leading T l ~ e term using an "NRG*" procedure involving 
NRG calculations for the spin-boson model supplemented 
by an ad hoc Hamiltonian term containing a coefficient 
that was adjusted until the inverse local susceptibility at 
the QCP had a temperature exponent close to \ . Figure 
[TOl provides a schematic comparison between the tem- 
perature dependences of the critical static local suscepti- 
bilities produced by the NRG (our data from Sec. MI Al 
above) and by the NRG* procedure (Figs. 7 and 8 of 
Ref. fish - If the leading term in x~H^) corresponded 
to x = ^ j then corrections arising from any subleading 
term with an exponent z > ^ would induce a crossover 
in x (T) around some T* to a higher-temperature be- 
havior with a greater slope. Instead, the e = 0.6 and 
0.7 NRG* results of Ref. [l8| appear to show a T z depen- 
dence with z < | for T > T*. One should expect, based 
on criticality, any such T z term to persist to sufficiently 
low temperatures that it, not the term, dominates 
the asymptotic behavior. 

The T X2 subleading term in the temperature depen- 
dence of the critical inverse static susceptibility is ac- 
companied by a subleading term in the frequency depen- 
dence of the inverse dynamical susceptibility of the form 
(—ioj) V2 . Our results are compatible with 1/2 = £2 and, in 
turn, with uj/T scaling for the subleading term (as well 
as the leading term, where y = x). If one interprets both 
the T x and T X2 terms in the temperature dependence as 
artifacts of NRG, then any intrinsic term must vary with 
a power of T even greater than X2', since x% > |, such 
behavior would be inconsistent with the validity of the 
quantum-to-classical mapping. Instead, the result rein- 
forces the conclusion that the fixed point is interacting 
and features uj/T scaling. 



8 



We now briefly discuss other critical exponents describ- 
ing the variation of the local magnetization M = (S z ), 
the corresponding susceptibility x — dM/dh\ h=0 , and 
the (imaginary-time) correlation length £ T in the vicinity 
of the QCP: 

M(T = 0, h = 0,g ></ c )oc (g - g c f , (18) 

M(T = 0,h,g = g c ) oc {h^ 6 , (19) 

X (T = 0,g<g c ) cx (g c -g)-\ (20) 

&K\g-g c \- v . (21) 

An important conclusion from the NRG treatment of the 
spin-boson model^ (as well as that of the Bose-Fermi 
Kondo model^) is that, together with x and y, these expo- 
nents satisfy hyperscalin g re lations 2 ^ expected to hold at 
an interacting QCP. Ref. Il8l argued that the NRG results 
for 8 and (3 are invalid for a reason that is entirely sep- 
arate from mass-flow error, namely the NRG truncation 
of the bosonic Hilbert space. (A variational procedure 
has recently been propose d 22 ' 23 to circumvent the effects 
of bosonic truncation.) However, it appears unnatural 
that two independent sources of errors should conspire 
to yield results that are consistent with hyperscaling. 



B. Path integral and Berry phase 

Any proper path integral requires taking the contin- 
uum limit in time. The quantum critical properties of 
the Bose-Fermi Kondo model or the spin-boson model 
can be described within a proper path integral repre- 
sentation. It was shown in Ref. H that the hyperscaling 
observed in the large-TV limit of the SU(7V) x SU(M) 
Bose-Fermi Kondo model even for ^ < e < 1 can be un- 
derstood as a consequence of a topological Berry-phase 
term that characterizes the path integral over SU(iV). 
Thus, the breakdown of the quantum-to-classical map- 
ping is not merely an artifact of the saddle point approx- 
imation, and it survives the inclusion of 1 /N corrections 4 
The path integral over the group SU(2) also involves a 
Berry-phase term that spoils the reinterpretation of the 
resulting complex action in terms of a classical action. In 
other words, the quantum-to-classical mapping cannot be 
applied for N ~ 2 and the quantum critical properties 
have to be obtained directly from the quantum model. 
A proper path integral on SU(2) is given by the func- 
tional integral on the Bloch sphere S 2 , the coset space of 
SU(2). Every path on S 2 is characterized by a geometric 
phase independent of the spin Hamiltonian. It therefore 
is natural to expect similar spin path integral represen- 
tations for the SU(2) Bose-Fermi Kondo model and for 
the spin-boson or easy-axis Bose-Fermi Kondo model. 

It was observed in Ref. that the construction of a 
Feynman path integral for the spin-boson model through 
a time-slicing procedure encounters difficulties in taking 
the continuum limit when using the orthonormal eigen- 
function basis of S z . The standard procedure for circum- 
venting this issue in the case of pure spin Hamiltonians 



is to rewrite the short-time propagator as a transfer ma- 
trix in spin spaced A generalization of this method to 
the spin-boson model rewrites the matrix element of the 
infinitesimal (imaginary) time evolution operator as 



|e-" T °|<r i+ K 



H+l) 



i+i J> 



(22) 



where _P l > I+1 is a matrix in the spin subspace, <f> is a c- 
number that labels the bosonic coherent states, e* > » < ^' i + 1 is 
related to the bosonic Berry phase, and i indexes the time 
slices. The matrix elements of P l ' l+1 are then expressed 
as the exponential of some function F(Si, <Si+i) with Si = 

zbltlft 



a a+b(Si+S, + 1 )+cSiS i+1 



Si,S i+1 



(23) 



with 



a = \ [ _CJT o ^t+i^i + m ( rT o)] , 

b= -\g-ro{<i>i+i + <i>i), 

c=\ [-ujtq 4>* i+1 <t>i - ln(rro)] . 



Inserting Eq. (|23|) into the expression for the partition 
function 



Z = Tre _/3ff 

M 

{<?M(pM\e- 

i=o o-i=t4 ' 

M-l 

x n 

k=Q 



M . 

—0 <T i= -t. ' 

M-l 



(24) 
Hn 'WoM 



(where A is a constant and t = j3/M), results in an ef- 
fective action having the form of the action for a classical 
spin chain coupled to the (/>-fields with nearest-neighbor 
interaction i [— utq 4>* + i4>i — ln(FTo)] , which is singular as 
To — > 0. However, the terms coupling the bosons are not 
singular for To — > 0, making it difficult to regularize this 
limit. Nonetheless, keeping track of both the transverse 
magnetic field Y and the coupling g to the bosonic bath 
is essential given that the competition between these two 
yields the QCP, and also that turning on g is a singular 
perturbation (as discussed at the end of Sec. IIV|) . A re- 
lated way to see this difficulty is that, for r — > 0, the 
Kondo-like energy scale goes to zero and the topological 
effect encoded in the spin flips is suppressed^ 

On the other hand, in a spin coherent-state represen- 
tation of the spin-boson model, the continuum limit fol- 
lowing a time-slicing procedure poses no difficulty. It 
leads to a well-defined path integral representation for 
the sub-ohmic spin-boson model [Eq. ([T])]^ 



Z = / V[n] exp[-Seff/2], 



(25) 
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where the effective action for the spin degrees of freedom 
is 



S e s = -iA[n] - / d,TYn x {T) 
Jo 

g 2 r? re 



(26) 



o Jo 



where 5 



+ *L / dT dT > nz{T)x -i {T _ T > )nz{T < )} 



cosh[w(/3/2 



Xo 1 (r-r)= / 



sinh(a;/3/2) 

(27) 

and n defines a point on S" 2 (|n(r)| 2 — 1). The effec- 
tive action of the sub-ohmic spin-boson problem in the 
continuum limit contains the Berry-phase term — i^4[n], 
where A[n] is equivalent to the area on the sphere S 2 
traced out by n(r) for < r < f3 with n(0) = n(/3). 
As discussed previously for models with continuous spin 
symmetry^, the Berry-phase term in the action, being 
imaginary, can invalidate the mapping of the quantum 
action to a classical one. 



VI. SUMMARY 

We have revisited the critical behavior of the sub- 
ohmic spin-boson problem. An analysis of the sublcading 
term in the temperature dependence of the inverse local 
spin susceptibility at the quantum critical point as calcu- 
lated using the numerical renormalization-group (NRG) 
method has provided evidence that the instrinsic behav- 
ior of the critical static spin susceptibility is given by 



Eq. [4] with a non-mean-field exponent. More specifically, 
we have examined the implications of our data for the 
two interpretations of the difference (for \ < e < 1) 
between Monte-Carlo calculations for a classical Ising 
spin-chain and NRG calculations for the quantum spin- 
boson model. The leading and subleading temperature 
dependences of the critical local susceptibility are consis- 
tent only with the interpretation that the classical spin- 
chain model fails to capture the critical behavior of the 
quantum model and that the criticality is described by 
an interacting fixed point characterized by lo/T scaling. 
This conclusion is further supported by the fact that the 
subleading frequency dependence of the inverse critical 
dynamical spin susceptibility has (within numerical un- 
certainty) the same exponent as that for the sublcading 
temperature dependence of the inverse critical static spin 
susceptibility. 



These results provide evidence for the violation of 
the quantum-to-classical mapping in the sub-ohmic spin- 
boson model. We have discussed the effect of a Berry- 
phase term in a continuum path-integral representation 
for this model. The importance of the Berry-phase term 
connects the violation of the quantum-to-classical map- 
ping observed here with that seen in models with contin- 
uous spin symmetry. 
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No. DMR-0710540, DMR-1006985, and DMR-1 107814, 
and by Robert A. Welch Foundation Grant No. C-1411. 
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